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Abstract 

We investigate how the addition of the singlet Higgs field in the NMSSM changes the 
sfermion branching ratios as compared to the MSSM. We concentrate in particular on the 
third generation, discussing decays of the heavier stop, sbottom or stau into the lighter 
mass eigenstate plus a scalar or pseudoscalar singlet Higgs. We also analyse stop, sbottom 
and stau decays into singlinos. It turns out that the branching ratios of these decays can 
be large, markedly influencing the sfermion phenomenology in the NMSSM. Moreover, we 
consider decays of first and second generation sfermions into singlinos. 

1 Introduction 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) provides an elegant solution 
to the fi problem of the MSSM by the addition of a gauge singlet superfield S [1]. The su- 
perpotential of the Higgs sector then has the form XS(Hd • H u ) + ^kS 3 . When S acquires a 
vacuum expectation value, this creates an effective /x term, /i = X(S), which is automatically 
of the right size, i.e. of the order of the electroweak scale. In this way, in the NMSSM the 
electroweak scale originates entirely from the SUSY breaking scale. 

The addition of the singlet field leads to a larger particle spectrum than in the MSSM: in 
addition to the MSSM fields, the NMSSM contains two extra neutral (singlet) Higgs fields -one 
scalar and one pseudoscalar- as well as an extra neutralino, the singlino. Owing to these extra 
states, the phenomenology of the NMSSM can be significantly different from the MSSM; not 
at least because the usual LEP limits do not apply to the singlet and singlino states. 

The NMSSM has recently become very popular. Most of the Feynman rules are given in [2]. 
The NMSSM Higgs phenomenology has been investigated extensively in [3-14] (a model variant 
without a S 3 term is discussed in [15]). Detailed studies of the neutralino sector are available 
in [16-23]. The relic density of (singlino) dark matter has been studied in [24]. The sfermion 
sector, on the other hand, has so far received very little attention, although here, too, one may 
observe differences as compared to the MSSM. In this letter we therefore investigate the decays 
of squarks and sleptons in the framework of the NMSSM and contrast them to the MSSM 
case. We concentrate on the third generation (stops, sbottoms and staus) where we expect the 
largest effects, but also consider decays of selectrons, smuons and lst/2nd generation squarks. 

In the MSSM, squarks and sleptons can decay via f i — > f i — > f'xf with i,j — 1, 2 (or 
i = L, R in case of no mixing) and k — 1, ...,4. Squarks can also decay into gluinos, q t — > qg, 



if the gluino is light enough. In addition, sfermions of the third generation (/ = t, b, r) can 
have the bosonic decay modes ~f i -> ft + W± and / 2 -»• f x + Z°, h°, H°, A [25]. In the 
NMSSM, we have additional decay modes into singlinos and singlet Higgs states: f i — > /x° 
with n = 1,...,5 and / 2 — > f x + A®, A®, H x , H®, H®. Pure singlets and singlinos couple in 
general very weakly to the rest of the spectrum. There are hence two potentially interesting 
cases: a ) large mixing of singlet and doublet Higgs states and/or large mixing of singlinos with 
gauginos-higgsinos and b) (very) light singlet /singlino states. In case b) A\ and H® are almost 
pure singlets with masses well below the LEP bound of rrih > 114 GeV, and the singlino is the 
LSP. We investigate these cases in this letter and show that they can markedly influence the 
sfermion phenomenology. 

The paper is organized as follows. In Section 2 we explain our notation, the potential and 
the relevant Feynman rules. In Section 3 we perform a numerical analysis, and in Section 4 we 
present our conclusions. 



2 Notation and couplings 
2.1 Potential 

We follow the notation of NMHDECAY [26]. The superpotential is then given as * 

W = h t Q- H u f c R + h b H d ■ QB C R + h T H d -Lr c R - XS H d ■ H u + ^kS 3 (1) 
with the SU (2) doublet superfields 




and the product of two SU(2) doublets 

X\ • X2 = X\X^ — X^X^ ■ (3) 
From Eq. (|1J) we derive the F-terms 



F Tl - 


htT R Hl - h h B c R H d , 


F'jpc — 


h t (T L H° u -B L H+), 


(4) 


Fb l = 


-h t T c R H+ + h b B R H° d , 


F B c = 


-h b (T L H d - B L H d ) , 


(5) 


Fr L = 


Kr c R H° d , 


F T c = 


-h T (v L H d -T L H° d ), 


(6) 


Fho = 


htTiT R — \SH d , 




—h t B L T R + \SH d , 


(7) 


F H o = 


h b B L B R + h T T L T R -\SH° u , 




-h b T L B' R -h T v L T c R + \SH+, 


(8) 



yielding the Yukawa part of the scalar potential = ^j^F^F*. Note that the F-terms in 
Eqs. flJJ) and (JSJ) imply direct interactions between the S field and the sfermions: 

V F D -h t X* (B L f c R S*H d * + f L f c R S*H° d *) - h\* (b l B c r S*H q u * + T L B C R S*H+^ 

h T \* (f L f c R S*H° u * + vjT R S*H+*) + h.c. . (9) 

1 Note the different signs of the h b and h T terms as compared to Eq. (A.l) of Ref. [26]. 

2 The superpotential Eq. Q posesses a discrete Z3 symmetry which is spontaneously broken at the electroweak 
phase transition. This results in cosmologically dangerous domain walls [27] . We implicitly assume a solution [28] 
to this domain wall problem which does not impact collider phenomenology. 
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We also need the soft SUSY-breaking potential for the derivation of the couplings, c.f. Eq. (A. 4) 
of [26], 

V sof t = h t A t Q ■ H u f c R + h b A b H d ■ QB C R + h T A T H d ■ Lr R - \A X S H d ■ H u + \kA k S 3 . (10) 
2.2 Sfermion— Higgs interaction 

In the following, we denote the neutral scalar and pseudoscalar Higgs bosons by Hf (i = 1,2,3) 
and A® (/ = 1,2), respectively. The interaction of Hf and A® with a pair of sfermions 
(j, k = 1,2) can be written as: 

^SH = dijkHffjfk+gfjk^-Ujfk- (11) 

Apart from D-term contributions, the Higgs-sfermion couplings are proportional to the Yukawa 
couplings hf. We therefore write g^ k explicitly for the third generation. For stops, we have 



(9ijkY - i^J^ ht ^*eff R )l R k2 + Mf R )2 Rt kl) - VdDjk^j S i2 

±= h t (A t R^Rl 2 + AIR^R^) - v u D jk + V2v u h 2 t 5 jk ^J S a 

+ -^j= v dh t (\*r]ir\ 2 + ^RtfRki) > 

(gf jk Y = -^hti^P^ + AtPn+v^P^R^Ri, 



(12) 



V2' 

+ ^=h t { Mf P l2 + A* t P n + vAP^Rl , (13) 
where /i e // is the effective \i term: 

fi eff = As (14) 

with s = (S) the vev of the singlet S. (In the presence of an additional generic \i term fj,H d H u , 
/i e ff — * n e ff = As + /!.) For sbottoms, we get 

& = ~ (j= h {A b R)lR b k2 + A* b R];Rl) + v d D jk + V2v d h 2 b 5 jk ^ S l2 

+ -±= v u h b (X*R i jl *Rl 2 + XR)* 2 R b kl ) S i3 , (15) 

(& = ~^=h b (A b P l2 + ti ff Pn + v u \*P lz )Ry i R b k2 

+ ^=h (A* b P l2 + M fPn + v u XP^)R)* 2 R b kl , (16) 

and analogously for staus with the obvious replacements h b — > h T , A b — > A T and R b — > R T . 
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T 3L 


Qt 


Yl 


Y R 


t 


1/2 


2/3 


1/3 


4/3 


b 


-1/2 


-1/3 


1/3 


-2/3 


f 


-1/2 


-1 


-1 


-2 



Table 1: Isospin, electric charge and hypercharges of stops, sbottoms and staus. 



In Eqs. (16)-(19), and are the Higgs mixing matrices as in [26], and W are the 
sfermion mixing matrices diagonalizing the sfermion mass matrices in the notation of [29,30]: 



h 
fx 



diag(mjk 



m 



2 

JL 

a f m f 



ajmf 

JR 



(17) 



where a t = A t — cot ft and a bjT = A h r — /z*^tan/3. Furthermore, = (H%), v u = (H®) and 



tan/3 = v u /vd- The D-terms are given by 



— ^ {{g 2 T 3 L - Y L g' 2 )R f ^R f kl - Y R g' 2 R f j2 * R f k2 



(18) 



where g and gr' are the SU(2) and £7(1) gauge couplings: g 2 = Ay/^GpM^ with Gf the Fermi 
constant and g' 2 = g 2 tan 2 = 4v / 2Gp(Mf — M^) using the on-shell relation sin 2 6 W = 
(1 — M^/M|) for the Weinberg angle 9yy] T 3L and Y L ^, are the 3rd component of the isospin 
and the hypercharge of the left (right) sfermion, respectively. We have Y = 2(Qf—T 3 ) where Qf 
is the electric charge. For completeness, the sfermion quantum numbers are listed in Table ^ 
Notice that in the CP-conserving case the pseudoscalars only couple to f 1 f 2 combinations 
and hence gf n = g[ 22 = in Eqs. (fT3j) and (|I5jl: moreover gf 12 = gf 2 i and g[ 12 = —g[ 2 i- Notice 
also that in the CP-violating case, the scalar and pseudoscalar Higgs states will mix to mass 
eigenstates h\ 5 similar to the MSSM case. The couplings of the charged Higgs bosons to 
sfermions are the same as in the MSSM. 



2.3 Sfermion— neutralino interaction 

The sfermion interaction with neutralinos has the same form as in the MSSM, 

%x° = 9 f («/> B + blP L ) xl fi + hx. . 
The only difference is the addition of the singlino state S: 

N nl B + N n2 W + N n3 H d + N ni H u + N n5 S 



(19) 



(20) 



with n = 1...5, N the matrix diagonalizing the 5x5 neutralino mass matrix in the basis 



M 



N 



(M x 

M 2 

-m z s w Cf3 mzcwcp 

mzswsp -m z c w sp 





-m z s w cp m z s w S/3 \ 

m z c w Cf3 -m z c w S{3 

-fi e ff -Ivsp 

-fJ-eff —IvCp 

—lvS/3 —lvC/3 2KS J 



(21) 



with sw = sin #14/, cw = cos6\y, sp = sin/3, cp = cos/3, and 



N*M.nN' = diag(m^o, m^a, m^o, m^o, m^o) . 



(22) 



We can hence use the couplings a( n and b in as given in [29,30] with the neutralino index running 
from 1 to 5 instead from 1 to 4. It is worth noting that the couplings between sfermions and 
singlinos only occur via the neutralino mixing. This is in contrast to the Higgs couplings where 
additional terms originating from Eq. Q are present in the sfermion-singlet interaction. Note 
also that [26] and [24] use the basis (B, W, H u , H^, S) for the neutralino system. Therefore, the 
indices 3 and 4 of the neutralino mixing matrix N need to be interchanged when comparing 
[24, 26] and this paper. 

3 Numerical results 

We have implemented all 2-body sparticle decays of the NMSSM in the SPheno [31] package. 
The Higgs sector of the NMSSM is calculated with NMHDECAY [26] linked to SPheno. For our 
discussion of sfermion decays in the NMSSM, we choose the benchmark scenario 3 of [24] which 
is characterized by 



with Mi = (g 1 /g 2 ) 2 M 2 ~ 0.5M 2 by GUT relations as an illustrative example. This leads to a 
light Xi with a mass of 35 GeV which is to 87% a singlino. The x 2 weights 169 GeV and is 
dominantly a higgsino. Moreover, we have a light scalar Higgs with a mass of = 36 GeV 
and a light pseudoscalar with rriA 1 = 56 GeV, both being almost pure singlet states and thus 
evading the LEP bounds. H®, H®, and A\ are SU(2) doublet fields similar to h°, H° and A 
in the MSSM. The relic density in this scenario is Qh 2 = 0.1155 [24]. 

In this letter, we are interested in the decays f 2 — > f\H®, f 2 — > f%A\, and / 12 — > fx®, 
with f = t, b, f. In order to see the relevance of these decays, we perform a random scan 
over the parameters of the third generation, Mq 3 , Mjy 3 , Mq 3 , Mg 3 , Mg 3 , A t , A&, A T . The 
sfermion mass parameters are varied between 100-800 GeV, and the trilinear couplings in their 
whole possible range allowed by the absence by charge or colour breaking minima. We compute 
the mass spectrum and the branching ratios at each scan point, accepting only points which 
pass the experimental bounds from LEP (the bounds from the LEP Higgs searches are fully 
implemented in NMHDECAY [26]). Owing to radiative corrections, the mass of H® var ies between 
~ 100 GeV and 117 GeV in the scan. The effect on the other quantities, in particular rriH 1 , 
rriA 1 and Qh 2 , is negligible. 

As a sideremark we note that renormalization-group (RG) arguments can be used to set 
lower bounds on the SU (2) doublet sfermion masses. Requiring, for example, that remain 

positive all the way up to the GUT scale implies m; > 0.9M 2 for the first and second generation 
at the weak scale. The corresponding bounds for the third generation are much lower because 
the Yukawa couplings contribute to the RG running with opposite sign as the SU{2) gauge 
coupling. We refrain, however, from imposing any such RG-inspired constraint in our analysis 
for two reasons: firstly because it is our aim to discuss the weak-scale phenomenology in 
the most general way using just one illustrative benchmark scenario, and secondly because the 



A = 0.4, k = 0.028, tan = 3, /V/ = As = 180 GeV, 
A x = 580 GeV, A K = -60 GeV, M 2 = 660 GeV, 



(23) 
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actual scale of SUSY breaking is unknown and may well be much lower than Mqvt (see e.g. the 
NMSSM variants of the models presented in [32]). 

Let us now discuss the sfermion branching ratios. Figure ^ shows scatter plots of the 
branching ratios of 62 and decays into A®, H® and H® as function of the heavy sfermion 
mass, m~ b2 or m? 2 . As can be seen, decays into the singlet Higgs bosons H® or A\ can have 
sizable branching ratios provided the sfermions are relatively light, mg ~ 2 < 400 GeV. This 
feature can easily be understood from the partial width of the decay of a heavier sfermion into 
a lighter one plus a massless singlet. For sbottoms we have, for instance, 



which is suppressed by a factor {v u /mi 2 ) 2 for heavy sbottoms. The factor c is c = cos 2 1d~ h for 
the scalar and c = 1 for the pseudoscalar singlet. Analogous expressions hold for staus with 
b — > t and for stops with b — > t and v u — > v<i- Decays into the SU{2) doublet Higgs can also 
have large branching fractions, provided the splitting of the sfermion mass eigenstates is large 
enough. Here note that for the parameter choice Eq. (j23J), ^4? and H® decay predominantly into 
bb and may hence only be distinguished by the different bb invariant masses. The H® on the 
other hand, decays to about 60% into H®H® and hence into a 46 final state. Both signatures, 
the one from decay into a Hf or A\ leading to bb with small invariant mass as well as the 46's 
from the decay into H®, are distinct from the usual MSSM case. 

We next turn to sfermion decays into singlinos. Figure 121 shows scatter plots of the branching 
ratios of t\ decays into txl, bxf, and of b\ decays into 6x12- As expected, t\ — > bxf and 
b\ — > bx® {Xi an d X2 being mostly doublet higgsinos) are in general the dominant modes if 
kinematically allowed. Nevertheless, as can be seen from the upper two plots of Fig. EJ decays 
into the singlino LSP can have sizable branching fractions, even if other decay modes are open. 
The size of the branching ratio into the singlino is mainly governed by its admixture from the 
doublet higgsinos, i. e. by the size of the A parameter. Similar features appear also in the stau 
decays as illustrated in Fig. El The pattern of f± is quite similar to that of bi, with > 10% 
branching ratio into the singlino for < 200-250 GeV, depending on the L/R character of 
the stau. For the ?2, the branching ratio of the decay into the singlino is even more important. 
In fact it can be 10-50% over a large part of the parameter space even if the decay into X2 * s 
open. 

In this context note also the xf and X2 wm cascade further into the singlino LSP, see 
e.g. [17,18,22]. This is quite distinct from the MSSM, where the singlino is absent and all 
decay chains end in what in our case is the Obviously, this also affects cascade decays of 
squarks of the first and second generation (and likewise of gluinos), since in case of a singlino 
LSP there is one more step in the chain as compared to the MSSM. At the LHC, a singlino LSP 
can in fact lead to similar signatures as a gravitino or axino LSP. The presence of light scalar or 
pseudoscalar Higgs bosons in the decay chains may be a way to distinguish the NMSSM from 
other scenarios. If A is large enough the decays f i — > fxi with 7^ NLSP may also be used for 
discrimination, since the corresponding f i decays into gravitino or axino would not occur. We 
will discuss this in more detail in a forthcoming paper. 

Last but not least we consider the decays of sfermions of the first two generations. Owing 
to the small Yukawa couplings, decays into Higgs bosons are negligible in this case. Decays 
into singlinos can, however, be important. As an example, Fig. 0] shows the branching ratios 
of e L R (being the same as those of £l l r ) for the scenario of Eq. (J2l?j) . For mg — 180-370 GeV, 




(24) 
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Figure 1: Branching ratios of b 2 (left) and r 2 (right) decays into Higgs bosons A®,!!®,!!® for 
scenario 3 of [24], c.f. Eq. ()23|). as function of the mass of the decaying particle. 



the decays e L — > v e Xi an d e R — > ex\ clearly dominate, giving a 2-step cascade decay into the 
singlino LSP. Nevertheless even in this case the decays e L R — > ex? have sizable rates, being of 
the order of 10% for cr. The reason is that the relative importance of the decays into charginos 
and neutralinos is determined by the gaugino components of the these particles, and \i an d 
X23 are mainly higgsino-like in our scenario. Only when the decay into x&i which is mainly a 
bino, gets kinematically allowed the direct decays into the singlino become negligible. Also for 
the decays of up and down squarks into singlinos we find branching ratios of 0(1-10)%. The 
implications for collider phenomenology will be discussed in detail elsewhere. 



4 Conclusions 

We have discussed the decays of sfermions in the NMSSM. We have shown that for stops, 
sbottoms and staus, in addition to the decay modes already present in the MSSM, decays into 
light singlet Higgs bosons, / 2 — > f x + H®, as well as decays into a singlino LSP, f i — > fxl 
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Figure 2: Branching ratios of t\ decays into tx± and bxf (left) and of &i decays into bx\ 
(right) for the parameters of Eq. /l ' v:M 



and 



(i = 1, 2) with x° — "S", can be important. This is in particular the case for light sfermions. 
The presence of these decay modes modifies the signatures of stop, sbottom and stau events 
as compared to the MSSM. Also for first and second generation sfermions it turned out that 
the decays into a singlino LSP can be quite important. Even if other decay modes are open, 
fi ~ * fXi with Xi — S can have (9(10%) branching ratio. Moreover, decays (of any SUSY 
particle) into singlinos or singlet Higgs bosons may significantly influence cascade decays of 
squarks and gluinos at the LHC. In particular in case of a singlino LSP, there is one more 
possible step in the decay chain than in the MSSM. A singlino LSP in the NMSSM can in fact 
lead to similar signatures at the LHC as a gravitino or axino LSP in the MSSM. The decays 
discussed in this letter may help discriminating these scenarios. 
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Figure 4: Branching ratios of £l (left) and £r (right) decays for the parameters of Eq. ([23]). 
The full, dashed, dotted, dash-dotted, and dash-dash-dotted lines are for the decays into ex?, 
ex®, ex®, ex4, and vxi > respectively. The branching ratios for smuons [Il,r are the same as for 
the selectrons. 



References 

[1] P. Fayet, Nucl. Phys. B 90 (1975) 104; H. P. Nilles, M. Srednicki and D. Wyler, Phys. 
Lett. B 120, 346 (1983); J. M. Frere, D. R. T. Jones and S. Raby, Nucl. Phys. B 222, 11 
(1983); J. P. Derendinger and C. A. Savoy, Nucl. Phys. B 237, 307 (1984); A. I. Veselov, 
M. I. Vysotsky and K. A. Ter-Martirosian, Sov. Phys. JETP 63, 489 (1986); J. R. Ellis, 
J. F. Gunion, H. E. Haber, L. Roszkowski and F. Zwirner, Phys. Rev. D 39, 844 (1989); 
L. Durand and J. L. Lopez, Phys. Lett. B 217, 463 (1989); M. Drees, Int. J. Mod. Phys. 
A 4, 3635 (1989). 

[2] F. Franke and H. Fraas, Int. J. Mod. Phys. A 12 (1997) 479 | arXiv:hep-ph/9512366 . 

[3] U. Ellwanger, Phys. Lett. B 303 (1993) 271 |arXiv:hep-ph/9302224| . 

[4] T. Elliott, S. F. King and P. L. White, Phys. Lett. B 31 4 (1993) 56 [arXiv:hep-ph/9305282| ; 
Phys. Rev. D 49 (1994) 2435 larXiv:hep-ph/9308309| . 

[5] U. Ellwanger, M. Rausch de Traubenberg and C. A. Savoy, Z. Phys. C 67 (1995) 665 
|arXiv:hep-ph/9502206| ; Nucl. Phys. B 492 (1997) 21 |arXiv:hep-ph/9611251| . 



9 



[6] S. F. King and P. L. White, Phys. Rev. D 52 (1995) 4183 |arXiv:hep-ph7 9505326|. 
[7] U. Ellwanger and C. Hugonie, Eur. Phys. J. C 25 (2002) 297 ||arXiv:hep-ph/9909260 . 



U. Ellwanger, J. F. Gunion and C. Hugonie, arXiv:hep-ph/0111179 l in "Les Houches 2001, 
Physics at TeV colliders" . 

[9] D. J. Miller, R. Ne vzorov and P. M. Zerwas, Nucl. Phys. B 681 (2004) 3 
|arXiv:hep"-ph/0304049|. 



[10] U. Ellwanger, J. F. Gunion, C. Hugonie and S. Moretti, arXiv:hep-ph/0305109, contribu- 
tion to the "LHC/LC Study Group Report". 



[11] U. Ellwanger, J. F. Gunion, C. Hugonie and S. Moretti, arXiv:hep-ph/0401228, in "Les 
Houches 2003, Physics at TeV colliders". 



[12] D. J. Miller and S. Moretti, |arXiv:hep-p h/0403137, contribution to the "LHC/LC Study 
Group Report". 

[13] U. Ellwanger, J. F. Gunion and C. Hugonie, |arXiv:hep-p h/0503203. 
[14] U. Ellwanger and C. Hugonie, |arXiv:hep- ph/0504269, 

[15] C. Panagiotakopoulos and A. Pilaftsis, Phys. Rev. D 63 (2001) 055003 



arXiv:hep-ph/0008268|. 



[16] F. Franke, H. Fraas and A. Bartl, Phys. Lett. B 336 (1994) 415 |arXiv:hep-ph/9408217 . 



[17] F. Franke and H. Fraas, Z. Phys. C 72 (1996) 309 | arXiv:hep-ph/9511275| . 

[18] U. Ellwanger and C. Hugonie, Eur. Phys. J. C 5 (1998) 723 | |arXiv:hep-ph/9712300| . 

[19] S. Hesselbach, F. Franke and H. Fraas, Phys. Lett. B 492 (2000) 140 
|arXiv:hep-ph/00073"T0| . 

[20] S. Hesselbach, F. F ranke and H. Fraas, Eur. Phys. J. C 23 (2002) 149 
|arXiv:hep-ph/01070"80| . 

[21] F. Franke and S. Hesselbach, Phys. Lett. B 526 (2002) 370 | arXiv:hep-ph/0111285| . 

[22] S. Y. Choi, D. J. Miller and P. M. Zerwas, Nucl. Phys. B 711 (2005) 83 
|arXiv:hep-ph/04072"09| . 



[23] G. Moortgat-Pick, S. Hesselbach, F. Franke and H. Fraas, |arXiv:hep-ph/ 0502036 



[24] G. Belanger, F. Boudjema, C. Hugonie, A. Pukhov and A. Semenov, 



arXiv:hep-ph/0505142, 

[25] A. Bartl et al, Phys. L ett. B 435 (199 8) 118 |arXiv:hep-ph/9804265| ; Phys. 

Lett. B 460 (1999) 1 57 |arXiv:hep-ph/9904417| ; Eur. Phys. J. directC 2 (2000) 6 
|arXiv:hep-ph/0002115| . 

[26] U. Ellwanger, J. F. Gunion and C. Hugonie, |arXiv:hep-ph70 406215 , 



10 



[27] S. A. Abel, S. Sar kar and P. L. White, Nucl. Phys. B 454 (1995) 663 
|arXiv:hep-ph/9506359| . 

[28] J. McDonald, Nucl. Phys. B 530 (1998) 325 |arXiv:hep-ph/9709512|; C. P anagiotakopoulos 
and K. Tamvakis, Phys. Lett. B 446 (1999) 224 |arXiv:hep-ph/9809475| . 

[29] A. Bartl, S. Hesselbach, K. Hidaka, T. Kernreiter and W. Porod, Phys. Rev. D 70 (2004) 
035003 |arXiv:hep-ph/0311338| . 

[30] T. Gajdosik, R. M. Godbole and S. Kraml, JHEP 0409 (2004) 051 |arXiv:hep-ph/0405167 . 

[31] W. Porod, Comput. Phys. Commun. 153 (2003) 275 larXiv:hep-ph/0301101| . The NMSSM 
extension will be published soon. In the meantime a private version can be obtained from 
W.P 

[32] B. C. Allanach, D Grellscheid and F. Quevedo, JHEP 0205 (2002) 048 
|arXiv:hep-ph/0111057 . 



11 



